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Abstract. We analyze the electronic Green's functions in the superconducting ground state of 
the t-J model using Gutzwiller-projected wave functions, and compare them to the conventional 
BCS form. Some of the properties of the BCS state are preserved by the projection: the total 
spectral weight is continuous around the quasiparticle node and approximately constant along the 
Fermi surface. On the other hand, the overall spectral weight is reduced by the projection with 
a momentum-dependent renormalization, and the projection produces electron-hole asymmetry in 
renormalization of the electron and hole spectral weights. The latter asymmetry leads to the bending 
of the effective Fermi surface which we define as the locus of equal electron and hole spectral 
weight. 
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INTRODUCTION 

High temperature superconductivity (HTSC) is one of tlie most intriguing phenomena 
in modern solid state physics. Experimentally, HTSC is observed in layered cuprate 
compounds. The undoped cuprates are antiferromagnetically ordered insulators which 
develop the characteristic superconducting "dome" upon doping with charge carriers. 

HTSC is interesting not only for promising technological applications, but also from 
a theoretical point of view. The relevant ingredients for HTSC are believed to be the 
following. 

• Low dimensionality (2d). 

• Strong short-range repulsion between electrons. 

• Doped Mott insulator. 

Taking these 3 ingredients, HTSC is modeled in the tight-binding description by large-i7 
Hubbard models or t-J models on the square lattice: 

HtJ = -tY.PG clcja Pg + JY. -Sj-"^) (1) 

i'J) (iJ) 

4-14- . 

where n = c^Cfj, S = jCaGcya'Ca' ^ ^re the Pauli matrices. The Gutzwiller projec- 
tor Pq = Tli{l — ni^riii^) prevents electrons from occupying the same lattice site. 



Repeated indices are summed over. 



The non-perturbative nature of the t-J model makes it an outstanding problem to 
solve in dimensions larger than one. Analytical techniques (renormalized or slave-boson 
mean-field theories [1,2]) are very crude and numerical techniques (e.g. exact diagonal- 
ization or cluster DMFT) are restricted to very small clusters or infinite dimensions, or 
they fail on the sign problem (QMC). An alternative approach was suggested by Ander- 
son shortly after the experimental discovery of HTSC, when he proposed a Gutzwiller- 
projected BCS wave function as superconducting ground state for cuprates [3]. Follow- 
ing this conjecture, many variational studies have been performed on the basis of what 
is called Anderson's (long range) RVB state. This state turned out to have very low vari- 
ational energy, close to exact ground state energies, as well as high overlap with the true 
ground states of small t-J clusters [4, 5]. On the other hand, many experimental facts 
about cuprate superconductors can be reproduced and are consistent with the variational 
results: e.g. clearly favored J-wave pairing symmetry, doping dependency of the nodal 
Fermi velocity and the nodal quasiparticle weight. Many of these successful efforts fol- 
lowing Anderson's proposal are summarized in the "plain vanilla RVB theory" of HTSC, 
recently reviewed in [6] . 

With help of the relatively recent technique of angle-resolved photoemission spec- 
troscopy (ARPES), experimentalists can probe the electronic structure of low-lying ex- 
citations inside the copper planes. The intensity measured in ARPES is proportional to 
the one-particle electronic spectral function: IpEs{k,0)) ocA{k,0)) [7]. It is therefore 
interesting to explore spectral properties within the framework of Gutz wilier-projected 
variational quasiparticle (QP) excitations. 

In this contribution we will discuss some of our results reported in [8]. For more 
details, in particular for more reference to experimental studies, we invite the reader to 
consult that paper. 

COHERENT SPECTRAL WEIGHTS 

Anderson's RVB state is given by 

\H)ocP„PQ\dBCS{A,fi)) . (2) 

We further define projected BCS quasiparticle excitations in a similar way, 

\H,k,o)<^PHPGrka\'^BCS) . (3) 

The unprojected states in Eqs. (2) and (3) are the usual ingredients of the BCS 
theory, |JBCS) = Uk^cYkalO), Jka = UtCta + (^nc^ks^ 4 = 5 " :|) = ^ " ^J' 

Ek = ^^^ + Al, = -2[cos{kj,) + cos{ky)] - [1, = A[cos{kj,) - cos (ky)]. Pq is the 

Gutzwiller projector and Ph projects on the subspace with H holes. \H) and \H, k, o) are 
normalized to one. The wave functions have two variational parameters, A and jU, which 
we adjusted to minimize the energy of the t-J Hamiltonian (1) for the experimentally 
relevant value 7 = 0.3 and every doping level. Note that in the RVB theory, A and 
jU are variational parameters without direct physical significance; physical quantities 



like excitation gap, superconducting order, or chemical potential must be calculated 
explicitly. 

The spectral weights of the coherent low-lying quasiparticles (3) can be written as 

Z+= \{H-l,k,a\cljH)\^ (4a) 
Z^-= \{H+l,k,G\c_i„am^- (4b) 

These weights are measured in ARPES experiments as the residues of the spectral 
function A (fc, (o) [7]. Note that in conventional BCS-theory, = and = v^. 

METHOD: VMC 

The variational Monte Carlo technique (VMC) allows to evaluate fermionic expectation 
values of the form for a given state In order to calculate the spectral 

weights (4) by VMC, the following exact relations can be used. 

Z+Z^- = \{H+l\ck^c_ki\H-l)W (5b) 

where x is hole doping [9, 8]. 

We use VMC to calculate the superconducting order parameter = I ('^fct'^-jfei) I' as 
well as diagonal matrix elements in the optimized t-J ground state (2). Using relations 
(5), we can then derive the spectral weights (4). The disadvantage of this procedure is 
large errorbars around the center of the Brillouin zone where both and 4>jfc are small. 

RESULTS 

In Fig. 1, we plot the nearest-neighbor superconducting orderparameter (c,|Cy|) as a 
function of doping. The curve shows close quantitative agreement with the result of 
Ref. [10], where the authors extracted the same quantity from the long-range asymp- 
totics of the nearest-neighbor pairing correlator, lim^^oo(cocgcjc|^^g). With the method 
employed here, we find the same qualitative and quantitative conclusions of previous 
authors [4, 10]: vanishing of superconductivity at half filling, x ^ 0, and at the su- 
perconducting transition on the overdoped side, Xc ^ 0.3. The optimal doping is near 
Xopt ~ 0.18. In the same plot we also show the commonly used Gutzwiller approxi- 
mation where the BCS orderparameter is renormalized by the factor gt = [1]. The 
Gutzwiller approximation underestimates the exact value by approximately 25%. 

In Fig. 2, we plot the spectral weights Z^, Z^, and Z^"' along the contour — > 
(0, ;r) ^ (;r, ;r) ^ in the Brillouin zone for different doping levels. Figure 3 shows 
the contour plots of Z^"' in the region of the Brillouin zone where our method produces 
small errorbars. From these data, we can make the following observations. 

• In the case of an unprojected BCS wave function, the total spectral weight is 
constant and unity over the Brillouin zone. Introducing the projection operator. 




X 

FIGURE 1. Doping dependency of the nearest-neighbor superconducting order parameter <I>,y (calcu- 
lated in the 14 x 14 system). The errorbars are smaller than the symbol size. The same quantity calculated 
in the Gutzwiller approximation is also shown for comparison. The variational parameter A is shown with 
the scale on the right. 



we see that for low doping (x ~ 3%), the spectral weight is reduced by a factor up 
to 20. The renormalization is asymmetric in the sense that the electronic spectral 
weight is more reduced than the hole spectral weight Z^. For higher doping 
(x ~ 23%), the spectral weight reduction is much smaller and the electron-hole 
asymmetry decreases. 
• Since there is no electron-hole mixing along the zone diagonal (J-wave), the spec- 
tral weights and Z^ have a discontinuity at the nodal point. Our data shows that 
the total spectral weight is continuous across the nodal point. Strong correlation 
does not affect these features of uncorrelated BCS-theory. 



Effective Fermi surface 

In strongly interacting Fermi systems, the notion of a Fermi surface (FS) is not clear 
at all. There are, however, several experimental definitions of the FS. Most commonly, 
kp is determined in ARPES experiments as the maximum of IVj^nj^l or the locus of 
minimal gap along some cut in the ^-plane. The theoretically better defined locus of 
nj^= 1/2 is also sometimes used. The various definitions of the FS usually agree within 
the experimental uncertainties [7]. The different definitions of the FS in HTSC were 
recently analyzed theoretically in Refs. [11, 12]. 

Here, we propose an alternative definition of the Fermi surface based on the ground 
state equal-time Green's functions. In the unprojected BCS state, the underlying FS is 
determined by the condition \uk\^ = |v^p. We will refer to this as the unprojected FS. 
Since \uk\^ and |vj;p are the residues of the QP poles in the BCS theory, it is natural to 
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FIGURE 2. QP spectral weights for 6 holes (upper left plot, x ~ 3%), 22 holes (lower left plot, jc ~ 11 %), 
34 holes (upper right plot, x ~ 17%), and 46 holes (lower right plot, x ~ 23%) on 196 sites. The spectral 
weights are plotted along the contour (0, n) {k, n) ^ Q (shown in inset). Plus signs (+) denote 
the spectral weight Z^, crosses (x) denote Z^, errorbars are shown. Solid dots denote their sum, the total 
spectral weight Z'^" , errorbars not shown. On the horizontal axis, the star (*) denotes the intersection with 
the unprojected Fermi surface along the — ^ (0, n) direction; the thick dot is the nodal point. Both Z^ 
and jump at the nodal point, while Z'^'' is continuous. The intersection with the effective Fermi surface 
(see text) is marked by an arrowhead. On the diagonal (last segment), k/ n is given in units of ^/2. 
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FIGURE 3. Contour plots of the total QP spectral weight Z'f. The effective FS (full line) and unpro- 
jected FS (dashed line) are also shown. The doping levels are x ~ 3%, 11%, 17%, and 23% (from left 
to right). The + signs indicate points where the values are known within small errorbars (see Section 
Method: VMC). 



replace them in the interacting case by and , respectively. We will therefore define 
the effective FS as the locus =Z^. 

In Fig. 3, we plot the unprojected and the effective FS which we obtained from VMC 
calculations. The contour plot of the total QP weight is also shown. It is interesting to 
note the following points. 

• In the underdoped region, the effective FS is open and bent outwards (hole-like 
FS). In the overdoped region, the effective FS closes and embraces more and more 



the unprojected one as doping is increased (electron-like FS). 

• Luttinger's rule [13] is clearly violated in the underdoped region, i.e. the area 
enclosed by the effective FS is not conserved by the interaction; it is larger than 
that of the unprojected FS. 

• In the optimally doped and overdoped region, the total spectral weight is approxi- 
mately constant along the effective FS within errorbars. In the highly underdoped 
region, we observe a small concentration of the spectral weight around the nodal 
point (~ 20%). 

Large "hole-like" FS in underdoped cuprates has also been reported in ARPES experi- 
ments by several groups [14, 15, 16]. 

It should be noted that a negative next-nearest hopping t' would lead to a similar 
FS curvature as we find in the underdoped region. We would like to emphasize that our 
original t-J Hamiltonian as well as the variational states do not contain any . Our results 
show that the outward curvature of the FS is due to strong repulsion, without need of t'. 
The next-nearest hopping terms in the microscopic description of the cuprates may not 
be necessary to explain the FS topology found in ARPES experiments. Remarkably, if 
the next-nearest hopping t' is included in the variational ansatz (and not in the original 
t-J Hamiltonian), a finite and negative t' is generated, as it was shown in Ref. [17]. 
Apparently, in this case the unprojected FS has the tendency to adjust to the effective FS. 
A similar bending of the FS was also reported in the recent analysis of the current carried 
by Gutz wilier-projected QPs [18]. A high-temperature expansion of the momentum 
distribution function of the t-J model was done in Ref. [19] where the authors find 
a violation of Luttinger's rule and a negative curvature of the FS. Our findings provide 
further evidence in this direction. 

A natural question is the role of superconductivity in the unconventional bending of 
the FS. In the limit A — > 0, the variational states are Gutz wilier-projected excitations 
of the Fermi sea and the spectral weights are step-functions at the (unprojected) FS. 
In a recent paper [20], it was shown that limj.^j^+Z^ = lim^^^^-Z^ for the projected 
Fermi-sea, which means that the unprojected and the effective FS coincide in that case. 
This suggests that the "hole-like" FS results from a non-trivial interplay between strong 
correlation and superconductivity. At the moment, we lack a qualitative explanation of 
this effect. However, it may be a consequence of the proximity of the system to the 
non- superconducting "staggered- flux" state [21, 22] or to antiferromagnetism near half- 
filling [10, 23]. 
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